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1. Elements of baryogenesis 

Particle physics unravels the structure of matter at short distances. The 
characteristic energy scales of the different layers of structure are of order 
1 eV, the binding energy of atoms, 1 MeV, the binding energy of nuclei, and 
100 GeV, the typical collision energy of quarks and leptons in present day 
high-energy colliders. Knowledge of the laws of nature which govern the 
interactions of elementary particles at these energies allows us to calculate 
the properties of a plasma of particles at the corresponding temperatures. 
Extrapolating the observed Hubble expansion of the universe back to early 
times one concludes that such temperatures must have been realized in the 
early stages of the evolution of the universe. 

At a temperature T ~ 1 eV electrons and nuclei combined to neutral 
atoms and the universe became transparent to photons. The observation 
and recent detailed investigation of the corresponding cosmic microwave 
background [1] is the basis of early-universe cosmology. The second main 
success is the prediction of the abundances of the light elements, D, 3 He, 
4 He and 7 Li based on the properties of the lightest baryons, protons and 
neutrons. The light elements were formed at a temperature T ~ 1 MeV. 
Agreement between theory and observation is obtained for a certain range 
of the parameter n, the ratio of baryon density and photon density [2], 

rj = — = (1.5 - 6.3) x 10" 10 , (1) 
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where the present number density of photons is n 7 ~ 400/cm 3 . Since no 
significant amount of antimatter is observed in the universe, the baryon 
density is equal to the cosmological baryon asymmetry, rig ~ - rig. 

The formation of the cosmological baryon asymmetry, i.e., the origin 
of the observed matter density today, can be understood in terms of the 
properties of quarks and leptons, the building blocks of the standard model 
which describes their interactons at energies of order 100 GeV and beyond. 
There are three families of left-handed and right-handed up-type and down- 
type quarks with baryon number B = 1/3 and lepton number L = 0, and 
three families of left-handed neutrinos and left- and right-handed electron- 
type leptons with B = and L = 1, 

QLi = ( ^ ) ' URi ' dm ' lu = ( et* ) ' eRi ' i = 1 ' ' ' 3 ' ^ 

Some of their interactions are not invariant under charge conjugation (C) 
and the combined charge conjugation and parity (CP) transformation [3]. 

A matter-antimatter asymmetry can be dynamically generated in an ex- 
panding universe if the particle interactions and the cosmological evolution 
satisfy Sakharov's conditions [4], i.e., 

— baryon number violation 

— C and CP violation 

— deviation from thermal equilibrium . 

At present there are a number of viable scenarios for baryogenesis. They can 
be classified according to the different ways in which Sakharov's conditions 
are realized. Already in the standard model C and CP are not conserved. 
Also B and L are violated by instanton processes [5]. In grand unified 
theories B and L are broken by the interactions of heavy gauge bosons 
and leptoquarks. This is the basis of the classical GUT baryogenesis [6]. 
Analogously, the L violating decays of heavy Majorana neutrinos lead to 
leptogenesis [7]. In supersymmetric theories the existence of approximately 
flat directions in the scalar potential provides new possibilities. Coherent 
oscillations of scalar fields may then generate large asymmetries [8]. 

The crucial departure from thermal equilibrium can also be realized in 
several ways. One possibility is a sufficiently strong first-order electroweak 
phase transition [9] . In this case CP violating interactions of the standard 
model or its supersymmetric extension could in principle generate the ob- 
served baryon asymmetry. However, due to the rather large lower bound 
on the Higgs boson mass of about 115 GeV, which is imposed by the LEP 
experiments, this interesting possibility is now restricted to a very small 
range of parameters in the supersymmetric standard model. In the case of 
the Affleck-Dine scenario the baryon asymmetry is generated at the end 
of an inflationary period as a coherent effect of scalar fields which leads 
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Figure 1. One of the 12-fermicm processes which are in thermal equilibrium in the 
high-temperature phase of the standard model. 



to an asymmetry between quarks and antiquarks after reheating [10]. For 
the classical GUT baryogenesis and for leptogenesis the departure from 
thermal equilibrium is due to the deviation of the number density of the 
decaying heavy particles from the equilibrium number density. How strong 
this deviation from thermal equilibrium is depends on the lifetime of the 
decaying heavy particles and the cosmological evolution. Further scenarios 
for baryogenesis are described in [11]. 

The theory of baryogenesis involves non-perturbative aspects of quan- 
tum field theory and also non-equilibrium statistical field theory, in partic- 
ular the theory of phase transitions and kinetic theory. A crucial ingredient 
is also the connection between baryon number and lepton number in the 
high-temperature, symmetric phase of the standard model. Due to the chi- 
ral nature of the weak interactions B and L are not conserved. At zero 
temperature this has no observable effect due to the smallness of the weak 
coupling. However, as the temperature approaches the critical temperature 
Tew of the electroweak transition, B and L violating processes come into 
thermal equilibrium [12]. 

The rate of these processes is related to the free energy of sphaleron- 
type field configurations which carry topological charge. In the standard 
model they lead to an effective interaction of all left-handed quarks and 
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Figure 2. Sphaleron rate in Bodeker's effective theory, two lattice implementations of 
HTL effective theory, and pure lattice theory interpreted as HTL effective theory [14] . 



leptons [5] (cf. fig. 1), 



B +L = ]J (iLiQLiQLilLi) , 



which violates baryon and lepton number by three units, 

AB = AL = 3 . 



(3) 



(4) 



The evaluation of the sphaleron rate in the symmetric high-temperature 
phase is a complicated problem. A clear physical picture has been obtained 
in Bodeker's effective theory [13] according to which low-frequency gauge 
field fluctuations satisfy the equation of motion 

D x B = ctE - C (5) 

Here ( is Gaussian noise, i.e., a random vector field with variance 

(C^xK^x')) = 2aS ij S(t - t')5(x - x') , (6) 

and a is a non-abelian conductivity. The sphaleron rate can then be written 
as [14], 

r~ (14.0 ±0.3)- {a w Tf . (7) 
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Figure 3. Effective lepton number violating interaction. 

A comparison with two lattice simulations is shown in fig. 2. From this 
one derives that B and L violating processes are in thermal equilibrium for 
temperatures in the range 

T EW ~ 100 GeV < T < T SPH ~ 10 12 GeV . (8) 

Sphaleron processes have a profound effect on the generation of the 
cosmological baryon asymmetry, in particular in connection with lepton 
number violating interactions between lepton and Higgs fields, 

£al=2 = 7^ A? ^iP c l Ljf + h.c. . (9) 

Such an interaction arises in particular from the exchange of heavy Ma- 
jorana neutrinos (cf. fig. 3). In the Higgs phase of the standard model, 
where the Higgs field acquires a vacuum expectation value, it gives rise to 
Majorana masses of the light neutrinos v e , and v T . 

Eq. (4) suggests that any B + L asymmetry generated before the elec- 
troweak phase transition, i.e., at temperatures T > Tew, will be washed 
out. However, since only left-handed fields couple to sphalerons, a non-zero 
value of B + L can persist in the high-temperature, symmetric phase if 
there exists a non-vanishing B — L asymmetry. An analysis of the chemi- 
cal potentials of all particle species in the high-temperature phase yields a 
relation between the baryon asymmetry Yb = {ub — ng)/s, where s is the 
entropy density, and the corresponding B — L and L asymmetries Yb~l 
and Yb, respectively [15], 

Y B = a Y B -l = Y L . (10) 

The number a depends on the other processes which are in thermal equi- 
librium. If these are all standard model interactions one has a = 28/79. If 
instead of the Yukawa interactions of the right-handed electron the AL = 2 
interactions (9) are in equilibrium one finds a = — 2/3 [16]. 
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The interplay between the sphaleron processes (fig. 1) and the lepton 
number violating processes (fig. 3) leads to an intriguing relation between 
neutrino properties and the cosmological baryon asymmetry. The decay 
of heavy Majorana neutrinos can quantitatively account for the observed 
asymmetry. 

2. Heavy particle decays in a thermal bath 

Let us now consider the simplest possibility for a departure from thermal 
equilibrium, the decay of heavy, weakly interacting particles in a thermal 
bath. To be specific, we choose the heavy particle to be a Majorana neutrino 
N = N c which can decay into a lepton Higgs pair l<p and also into the CP 
conjugate state T<j> 



In the case of CP violating couplings a lepton asymmetry can be generated 
in the decays of the heavy neutrinos N which is then partially transformed 
into a baryon asymmetry by sphaleron processes [7]. Compared to other sce- 
narios of baryogenesis this leptogenesis mechanism has the advantage that, 
at least in principle, the resulting baryon asymmetry is entirely determined 
by neutrino properties. 



Figure 4- AL = 1 processes: decays and inverse decays of a heavy Majorana neutrino. 



N —> I </> , N —>l (/> . 



(11) 





Figure 5. AL = and AL = 2 lepton Higgs processes. 



2.1. BOLTZMANN EQUATIONS 

The generation of a baryon asymmetry is an out-of-equilibrium process 
which is generally treated by means of Boltzmann equations. A detailed 
discussion of the basic ideas and some of the subtleties has been given in 
[17]. The main processes in the thermal bath are the decays and the inverse 
decays of the heavy neutrinos (cf. fig. 4), and the lepton number conserving 
(AL = 0) and violating (AL = 2) processes (cf. fig. 5). In addition there 
are other processes, in particular those involving the t-quark, which are also 
important in a quantitative analysis [18, 19]. A lepton asymmetry can be 
dynamically generated in an expanding universe if the partial decay widths 
of the heavy neutrino violate CP invariance, i.e., 



Here T is the total decay width, and the parameter e ^ 1 measures the 
amount of CP violation. 

The Boltzmann equations for the number densities of heavy neutrinos 
(njv), leptons (ri/) and antileptons (nr) corresponding to the processes in 
figs. 4 and 5 are given by 




(12) 



driN 



+ 3HriN 



-j(N -> l<p) + y(l(f> -> N) 
- 7 (iV 14>) + 7(Z0 N) , 
7 (iV -» l4>)-i(l(f> -► N) 
+ 7 (ty - l<f>) ~ - U) , 

7 (iv -» 14>) - ~/(i4> N) 
+ 7 (14> - U) - 7 (/0 l<t>) , 



dt 



(13) 




(14) 




(15) 
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Figure 6. Time evolution of the number density to entropy density ratio. AtT^M the 
system gets out of equilibrium and an asymmetry is produced. 

with the reaction rates 

7 (jV-ty) = J d$ 123 f N { Pl )\M(N -^l<t>)\ 2 (16) 
7 (ty _> ty) = / d$ 1234 /i(pi)^(p2)|-M'(Z^ - /^)| 2 , • • • (17) 

Here H is the Hubble parameter, d$i... n denotes the phase space integration 
over particles in initial and final states, 

n = d3 f 1 . . . d3 f" (27r) 4 <5 4 ( Pl + . . . - p n ) , (18) 

(2tt) 3 2£i (27r) 3 2£ n v ; VF F7 v ' 

and 

/i(p) = exp {-PEiip)) , ^(p) = 9i J j^fiip) , (19) 

are Boltzmann distribution and number density of particle i = N,l,<p at 
temperature T = 1/(3, respectively. M. and M.' denote the scattering matrix 
elements of the indicated processes at zero temperature; the prime indicates 
that for the 2^2 processes the contribution of the intermediate resonance 
state has been subtracted. For simplicity we have used in eqs. (16) and 
(17) Boltzmann distributions rather than Bose-Einstein and Fermi-Dirac 
distributions, and we have also neglected the distribution functions in the 
final state, which is a good approximation for small number densities. Sub- 
tracting (15) from (14) yields the Boltzmann equation for the asymmetry 

A typical solution of the Boltzmann equations (13) - (15) is shown in 
fig. 6. Here the ratios of number densities and entropy density, 




s 



(20) 
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are plotted, which remain constant in an expanding universe in thermal 
equilibrium. A heavy neutrino, which is weakly coupled to the thermal 
bath, falls out of thermal equilibrium at temperatures T ~ M since its 
decay is too slow to follow the rapidly decreasing equilibrium distribution 
/at ~ exp(— (3M). This leads to an excess of the number density, > n$. 
CP violating partial decay widths then yield a lepton asymmetry which, 
by means of sphaleron processes, is partially transformed into a baryon 
asymmetry. 

The Boltzmann equations are classical equations for the time evolution 
of number densities. The collision terms, however, are 5-matrix elements 
which involve quantum mechanical interferences of different amplitudes in 
a crucial manner. Since these scattering matrix elements are evaluated at 
zero temperature, one may worry to what extent the quantum mechani- 
cal interferences are affected by interactions with the thermal bath. A full 
quantum mechanical treatment can be based either on the density matrix 
[20] or on on the Kadanoff-Baym equations [21]. 

Another subtlety is the separation of the 2^2 matrix elements into a 
resonance contribution and remainder [17], 

\M{l<p - U)\ 2 = \M'{l<p - U)\ 2 + \M res (l4> - W , (21) 

where the resonance contribution has the form 

M T es{l<t> U) oc M{l<j) N)M{N Uf = \M{l<j) A0| 2 • (22) 

The entire effect of baryon number generation crucially depends on this 
separation. The particles which participate in the 2^2 processes are 
massless, hence their distribution functions always coincide with the equi- 
librium distribution. Only the resonances, treated as on-shell particles, fall 
out of thermal equilibrium and can therefore generate an asymmetry in 
their decays. General theoretical arguments require cancellations between 
these two types of contributions which is illustrated by the following exam- 
ple. 

2.2. CANCELLATIONS IN THERMAL EQUILIBRIUM 

If all processes, including those which violate baryon number, are in thermal 
equilibrium the baryon asymmetry vanishes. This is a direct consequence 
of the CPT invariance of the theory, 

(B) = Tr(pB) = Tr((CPT){CPT)- 1 exp{-pH)B^) 

= Tr (exp {-pH^CPT^BiCPT)) = -Tv(pB) = . (23) 
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Hence, no asymmetry can be generated in equilibrium, and the transition 
rate which determines the change of the asymmetry has to vanish, 

d{ni ^ t ni) + 3H( ni - nj) = M eq = , (24) 

where the superscript eq denotes rates evaluated with equilibrium distribu- 
tions. 

From eqs. (12), (14) and (15) one obtains for the resonance contribution, 
i.e. decay and inverse decay, 

A 7r e f s = -2e^(N -+ 14>) . (25) 

This means in particular that the asymmetry generated in the decay is 
not compensated by the effect of inverse decays. On the contrary, both 
processes contribute the same amount. 

The rate A 7 ^f s has to be compensated by the contribution from 2^2 
processes which is given by 

A 72 e i 2 = 2 j d<s> 12M ft q (pi)f; g ( P 2) (ww - m\ 2 - \m'{T4> - m 2 ) . 

(26) 

For weakly coupled heavy neutrinos, i.e. T oc A 2 M with A 2 < 1, this 
compensation can be easily shown using the unitarity of the 5-matrix. 
The sum over states in the unitarity relation, 

]T ({MM - X)\ 2 - \M(X -+ lcf>)\ 2 ) = , (27) 
x 

can be restricted to two-particle states to leading order in the case of weak 
coupling A. This implies for the considered 2^2 processes, 

Yl (\M(l<f> - /W - \M(T'& - m 2 ) = , (28) 

l,cp,l',<f>' 

where the summation includes momentum integrations under the con- 
straint of fixed total momentum. From eqs. (21) and (28) one obtains 

A 72 e i 2 = 2 J ' d^ft q ( Pl )fl q {p2) 

(HAW^W + IAW^^)! 2 ) • (29) 

In the narrow width approximation, i.e. to leading order in A 2 , this yields 
the wanted result, 

= 2 J d^ 12U f^( Pl )r 4> \p 2 ) {-\M(l<f> - N)\ 2 \M{N - U)? 

+ \M(U - iV)| 2 |^(iV - m 2 ) ^6(8 - M 2 ) 
= 2e^(N -+ l<f>) = -A 7r e f s . ' (30) 
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This cancellation also illustrates that the Boltzmann equations treat res- 
onances as on-shell real particles. Off-shell effects can only be taken into 
account in a full quantum mechanical treatment. 



3. Cosmological bounds on neutrino masses 

Leptogenesis requires lepton number violation. On the other hand, the 
existence of a non-vanishing baryon asymmetry also restricts the allowed 
amount of lepton number violation and implies upper bounds on the Ma- 
jorana masses of the light neutrinos. 

In the standard model (SM) neutrinos are massless. In general, however, 
the exchange of heavy particles gives rise to an effective lepton Higgs in- 
teraction which, after electroweak symmetry breaking, generates neutrino 
masses. The lagrangian describing all fermion-Higgs couplings then reads 



£>Y = -h di jd,RiqLjHi - h uij URiqLjH 2 - h ei jeRilLjHi 
1 
2 



+^fij luH 2 C l Lj H 2 + h.c. . (31) 



Here qn, URi, dju, In, eRi, i = 1...N, are N generations of quark and 
lepton fields, H\ and H 2 are Higgs fields with vacuum expectation values 
V{ = (Hf) ^ 0. ha, h u , h e and / are NxN complex matrices. For the further 
discussion it is convenient to choose a basis where h u and h e are diagonal 



and real. In the case v\ ~ v 2 ~ v = y v\ + v\ the smallest Yukawa coupling 
is h e n = rn e /v\ followed by h u n = m u /v 2 . 

The mixing of the different quark generations is given by the Kobayashi- 
Maskawa matrix Vd, which is defined by 

hd Vd = — ■ (32) 

Vi 

Here rrid is the diagonal real down quark mass matrix, and the weak eigen- 
states of the right-handed d-quarks have been chosen to be identical to the 
mass eigenstates. Correspondingly, the mixing matrix V v in the leptonic 
charged current is determined by 

V? f V„ = - ^ , (33) 

where 

/ mi \ 

m u = m 2 (34) 

\ m 3 / 

is the diagonal and real mass matrix of the light Majorana neutrinos. 
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3.1. CHEMICAL EQUILIBRIUM 

In a weakly coupled plasma with temperature T and volume V one can 
assign a chemical potential fx to each of the quark, lepton and Higgs fields. 
In the SM with one Higgs doublet, i.e., H2 = Hi = if, and N generations 
one has 5N+1 chemical potentials. The corresponding partition function is 
[22] 

Z(fi, T, V) = Tre~^ H ~^i ^ Qi) . (35) 

Here (3 = 1/T, H is the Hamilton operator and Qi are the charge operators 
for the corresponding fields. The asymmetry in the particle and antiparticle 
number densities is then given by the derivative of the thermodynamic 
potential, 

n i -n i = - dn ^ T) , Sl(j i ,T) = ~hxZ( Ji ,T,V). (36) 
For a non-interacting gas of massless particles one has 

gT 3 \ Pm + O ((/3/ij) 3 ) , fermions , 

n; — rii = < > ,< (37) 

6 1 2Pm + 0((Pm) 3 ) , bosons. 



The following analysis will be based on these relations for (3fii <C 1. How- 
ever, one should keep in mind that the plasma of the early universe is 
very different from a weakly coupled relativistic gas due to the presence of 
unscreened non-abelian gauge interactions. Hence, non-perturbative effects 
may be important in some cases. 

In the high-temperature plasma quarks, leptons and Higgs bosons in- 
teract via Yukawa and gauge couplings and, in addition, via the non- 
perturbative sphaleron processes. In thermal equilibrium all these processes 
yield constraints between the various chemical potentials. The effective in- 
teraction (3) induced by the SU (2) electroweak instantons yields the con- 
straint [12], 

E ( 3 /V + A*k) = • (38) 

i 

One also has to take the SU(3) QCD instanton processes into account [23] 
which generate the effective interaction 

A = \{(q mLt u c m d c Ri ) (39) 

i 

between left-handed and right-handed quarks. The corresponding relation 
between the chemical potentials reads 

{fyqi ~ V-ui ~ Mdi) = . (40) 
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A third condition, which is valid at all temperatures, arises from the re- 
quirement that the total hypercharge of the plasma vanishes. From eq. (37) 
and the known hypercharges one obtains 

(ftqi + 2>Hui ~ Hdi ~ Mi - Mi + = " ( 41 ) 

i 

The Yukawa interactions, supplemented by gauge interactions, yield re- 
lations between the chemical potentials of left-handed and right-handed 
fermions, 

Mji _ /V ~ Vdj = , fl qi + Mp~ Vuj = , Mi~ Mp~ Mj = • (42) 

Furthermore, the AL = 2 interaction in (31) implies 

Mi + fj,<p = ■ (43) 

The above relations between chemical potentials hold if the correspond- 
ing interactions are in thermal equilibrium. In the temperature range Tew ~ 
100 GeV < T < T SPH ~ 10 12 GeV, which is of interest for baryogenesis, 
this is the case for all gauge interactions. It is not always true, however, 
for Yukawa interactions. The rate of a scattering process between left- and 
right-handed fermions, Higgs boson and W-boson, 

^LV iPrW , (44) 

is T ~ a\ 2 T, with a = g 2 /(4ir). This rate has to be compared with the 
Hubble rate, 

The equilibrium condition T(T) > H(T) is satisfied for sufficiently small 
temperatures, 

T < T x ~ A 2 10 16 GeV . (46) 

Hence, one obtains for the decoupling temperatures of right-handed elec- 
trons, up-quarks,..., 

T e ~ 10 4 GeV , T u ~ 10 6 GeV , . . . . (47) 

At a temperature T ~ 10 10 GeV, which is characteristic of leptogenesis, 
eR = em, hr = e R2 , d R = dm, SR = d R2 and u R = u R1 are out of 
equilibrium. 

Note that the QCD sphaleron constraint (40) is automatically satisfied if 
the quark Yukawa interactions are in equilibrium (cf. (42)). If the Yukawa 
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interaction of one of the right-handed quarks is too weak, the sphaleron 
constraint still establishes full chemical equilibrium. 

Using eq. (37) also the baryon number density rig = BT 2 /6 and the 
lepton number densities iil = LT 2 /Q can be expressed in terms of the 
chemical potentials. The baryon asymmetry B and the lepton asymmetries 
Li read 

B = Yl + Mu* + Vdi) , (48) 

i 

Li = 2fi H + n e i, L = ^Li. (49) 



3.2. RELATIONS BETWEEN B, L AND B — L 

Knowing which particle species are in thermal equilibrium one can derive 
relations between different asymmetries. Consider first the most familiar 
case where all Yukawa interactions are in equilibrium and the AL = 2 
lepton-Higgs interaction is out of equilibrium. In this case the asymmetries 
Li — B/N are conserved. The Yukawa interactions establish equilibrium 
between the different generations, 

Hi = W , Vqi = , etc. (50) 

Together with the sphaleron process and the hypercharge constraint they 
allow to express all chemical potentials, and therefore all asymmetries, in 
terms of a single chemical potential which may be chosen to be //j. The 
result reads 

2N + 3 6N + 1 _ 2N - 1 

^ ~ ^v + 3 W ' Md "~6iV + 3 W ' Mu ~6F+3 W ' 

Vq = -3W , = 6iv + 3 W • (51) 

The corresponding baryon and lepton asymmetries are 

4N UN 2 + 9N 
B = iii > L = — at , (52) 

which yields the well-known connection between the B and B — L asym- 
metries [24] 

B= 8N + A (B-L). (53) 
22N + 13 V ; v ; 

Note, that this relation only holds for temperatures T 3> v. In general, the 
ratio B/(B - L) is a function of v/T [25, 26]. 
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Another instructive example is the case where the AL = 2 interactions 
are in equilibrium but the right-handed electrons are not. Depending on the 
neutrino masses and mixings, this could be the case for temperatures above 
T e ~ 10 4 GeV [27]. Right-handed electron number would then be conserved, 
and Yukawa and gauge interactions would relate all asymmetries to the 
asymmetry of right-handed electrons. The various chemical potentials are 
given by (/x e = fi el , /j, s = fi e2 = . . . = /x e jv), 

3 1 1 

Ml = "^e, ^ = ^e, *P=^*>- (54) 

The corresponding baryon and lepton asymmetries are [27] 

1 4iV + 3 

B = -fi e , L= wN ii e , (55) 

which yields for the relation between B and B — L, 

2N 

B = (B-L). (56) 

2iV + 3 v ; v ' 

Note that although sphaleron processes and AL = 2 processes are in equi- 
librium, the asymmetries in B, L and B — L do not vanish! 



3.3. CONSTRAINT ON MAJORANA NEUTRINO MASSES 

The two examples illustrate the connection between lepton number and 
baryon number induced by sphaleron processes. They also show how this 
connection depends on other processes in the high-temperature plasma. 
To have one quark-Higgs or lepton-Higgs interaction out of equilibrium is 
sufficient in order to have non- vanishing B, L and B — L. If all interactions 
in (31) are in equilibrium, eqs. (43) and (51) together imply m = and 
therefore 

B = L = B - L = , (57) 

which is inconsistent with the existence of a matter-antimatter asymmetry. 
Since the equilibrium conditions of the various interactions are temperature 
dependent, and the AL = 2 interaction is related to neutrino masses and 
mixings, one obtains important constraints on neutrino properties from the 
existence of the cosmological baryon asymmetry. 

The AL = 2 processes described by (9) take place with the rate [28] 

r A L= 2 (T) = i^ E <■ (58) 

i=e,fi,r 
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Requiring T/^l = 2(T) < H(T) then yields an upper bound on Majorana 
neutrino masses, 

£<<(o.2eV (^) V2 ) 2 . (59) 

For typical leptogenesis temperatures T ~ 10 10 GeV this bound is com- 
parable to the upper bound on the electron neutrino mass obtained from 
neutrinoless double beta decay. Note, that the bound also applies to the 
r-neutrino mass. However, if one uses for T the decoupling temperature of 
right-handed electrons, T e ~ 10 4 GeV, the much weaker bound m u < 2 keV 
is obtained [27]. 

Clearly, what temperature one has to use in eq. (59) depends on the ther- 
mal history of the early universe. Some information is needed on what kind 
of asymmetries may have been generated as the temperature decreased. 
This, together with the temperature dependence of the lepton-Higgs inter- 
actions then yields constraints on neutrino masses. 

3.4. PRIMORDIAL ASYMMETRIES 

The possible generation of asymmetries can be systematically studied by 
listing all the higher-dimensional SU(3) x SU(2) x U(l) operators which 
may be generated by the exchange of heavy particles. The dynamics of the 
heavy particles may then generate an asymmetry in the quantum numbers 
carried by the massless fields which appear in the operator. 

For d=5, there is a unique operator, which has already been discussed 
above, 

(!u<p)(lLj<p) ■ (60) 

It is generated in particular by the exchange of heavy Majorana neutrinos 
whose coupling to the massless fields is 

KijVRihjV ■ (61) 

The out-of-equilibrium decays of the heavy neutrinos can generate a lep- 
ton asymmetry, which is the well-known mechanism of leptogenesis. The 
decays yield asymmetries Li — B/N which are conserved in the subsequent 
evolution. The initial asymmetry in right-handed electrons is zero. In order 
to satisfy the out-of-equilibrium condition it is very important that at least 
some Yukawa couplings are small and that the right-handed neutrinos carry 
no quantum numbers with respect to unbroken gauge symmetries. 

In order to study possible asymmetries of right-handed electrons one has 
to consider operators containing e#. A simple example, with d=6, reads 



(qLilLj){u c Rk e c Rl ) 



(62) 
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It can be generated by leptoquark exchange (% ~ (3*, 1, 1/3)), 

(\%q T Lt Cl Ll +\^u T m Ce R] ) X . (63) 

Note, that <p and % form a 5*-plet of .St/ (5). In principle, out-of-equilibrium 
decays of leptoquarks may generate a en asymmetry. One may worry, how- 
ever, whether the branching ratio into final states containing en is suf- 
ficiently large. Furthermore, it appears very difficult to satisfy the out- 
of-equilibrium condition since leptoquarks carry colour. Maybe, all these 
problems can be overcome by making use of coherent oscillations of scalar 
fields [8] or by special particle production mechanisms after inflation. How- 
ever, we are not aware of a consistent scenario for the generation of a 
ei?-asymmetry. Hence, it appears appropriate to take the bound eq. (59) 
on Majorana neutrino masses as a guideline and to examine its validity in 
each particular model. 

4. Neutrino masses and mixings 

Majorana masses for the light neutrinos are most easily generated by the 
exchange of heavy Majorana neutrinos. Such heavy 'right-handed' neutrinos 
are predicted by all extensions of the standard model which contain B — L 
as a local symmetry. The most general Lagrangian for couplings and masses 
of charged leptons and neutrinos reads 

C Y = -heijeRilLjHx - h l/ij VR i l Lj H 2 - ^h rij u c Ri u Rj R + h.c. . (64) 

The vacuum expectation values of the Higgs fields, (H\) = v\ and (H2) = 
V2 = tan/3 vi, generate Dirac masses m e and mo for charged leptons and 
neutrinos, m e = h e v\ and mo = h v V2, respectively, which are assumed to 
be much smaller than the Majorana masses M = h r (R). This yields light 
and heavy neutrino mass eigenstates according to the seesaw mechanism 
[29], 

v~Vjv L + v c L VZ , + (65) 

with masses ^ 

m v ~ -V^m^—m D V v , m N ~ M . (66) 

Here V v is the mixing matrix in the leptonic charged current (cf. eqs. (31)- 
(33)). 

In models of leptogenesis the predicted value of the baryon asymmetry 
depends on the CP asymmetry (cf. (12)) which is determined by the Dirac 
and the Majorana neutrino mass matrices. Depending on the neutrino mass 
hierarchy and the size of the mixing angles the CP asymmetry can vary over 
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many orders of magnitude. It is therefore important to see whether patterns 
of neutrino masses [30] motivated by other considerations are consistent 
with leptogenesis. In the following we shall consider two examples. 

An attractive framework to explain the observed mass hierarchies of 
quarks and charged leptons is the Froggatt-Nielsen mechanism [31] based 
on a spontaneously broken U(1)f generation symmetry. The Yukawa cou- 
plings arise from non-renormalizable interactions after a gauge singlet field 
<I> acquires a vacuum expectation value, 



Here gij are couplings 0(1) and Qi are the U(1)f charges of the various 
fermions, with Q$ = — 1. The interaction scale A is usually chosen to be 
very large, A > Aqut- In the following we shall discuss two different re- 
alizations of this idea which are motivated by the atmospheric neutrino 
anomaly [32]. Both scenarios have a large — u T mixing angle. They dif- 
fer, however, by the symmetry structure and by the size of the parameter 
e which characterizes the flavour mixing. 

4.1. SU(5) x U(1) F 

This symmetry has been considered by a number of authors. Particularly 
interesting is the case with a non-parallel family structure where the chiral 
U(1)f charges are different for the 5*-plets and the 10-plets of the same 
family [33]-[37]. An example of possible charges Qi is given in table 1. 

The assignment of the same charge to the lepton doublets of the second 
and third generation leads to a neutrino mass matrix of the form [33, 34], 



This structure immediately yields a large z/„ — v T mixing angle. The phe- 
nomenology of neutrino oscillations depends on the unspecified coefficients 




(67) 




(68) 



i>i e R3 e R2 e Rl lL3 lL2 I LI V°R3 ^R2 v RX 







1 



2 











1 







1 



2 



TABLE 1. Chiral charges of charged and neutral leptons with 
SU(5) x U(1) F symmetry [37]. 
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0(1). The parameter e which gives the flavour mixing is chosen to be 

which corresponds to the mass ratio m^/m T . The three Yukawa matrices 
for the leptons are given by, 

e 3 e 2 e 2 \ / e 4 e 3 e 2 

h e , K ~ | e 2 e e , /i r ~ e 3 e 2 e | . (70) 
ell/ V e 2 e 1 ' 

Note, that h e and /i^ have the same, non-symmetric structure. One easily 
verifies that the mass ratios for charged leptons, heavy and light Majorana 
neutrinos are given by 

m e : : m T ~ e 3 : e : 1 , M x : M 2 : M 3 ~ e 4 : e 2 : 1 , (71) 

mi : ni2 : ~ e 2 : 1 : 1 . (72) 

The masses of the two eigenstates v 2 and v-$ depend on the unspecified 
factors of order one, and may differ by an order of magnitude [38, 39]. They 
can therefore be consistent with the mass differences AmJ w ~ 4- 10~ 6 — 1 • 
10~ 5 eV 2 [40] inferred from the MSW solution of the solar neutrino problem 
[41] and Ara^ 3 ~ (5 • 10~ 4 — 6 • 10~ 3 ) eV 2 associated with the atmospheric 
neutrino deficit [32] . In the following we we shall use for numerical estimates 
the average of the neutrino masses of the second and third family, m u = 
)V2 „ 10-2 eV. 

The choice of the charges in table 1 corresponds to large Yukawa cou- 
plings of the third generation. For the mass of the heaviest Majorana neu- 
trino one finds 

v 2 

M 3 ~ ^- ~ 10 15 GeV . (73) 
m v 

Since /i r 33 and the gauge coupling of U(1)b-l are 0(1)> this implies that 
B — L is broken at the unification scale Aqut- 



4.2. SU(3) C x SU(3) L x SU{3) R x U(l) F 

This symmetry arises in unified theories based on the gauge group E^- 
The leptons e c R , II and v c R are contained in a single (1, 3, 3) representation. 
Hence, all leptons of the same generation have the same U(1)f charge and 
all leptonic Yukawa matrices are symmetric. Masses and mixings of quarks 
and charged leptons can be successfully described by using the charges given 
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i>i e R3 e R2 e Rl lL3 lh2 I LI Vrz V R2 v h 



Qi 



o 



o 







TABLE 2. Chiral charges of charged and neutral leptons with 
SU(3) C x SU(3)l x SU(3)r x U{1)f symmetry [36]. 



in table 2. Clearly, all three Yukawa matrices have the same structure 1 , 



h e , h r ~ 



V 



e 

F 5/2 



£ 3 £ 5/2 \ 

eV2 l y 



V 



e 



e 

cl/2 



1 / 



, (74) 



but the expansion parameter in /i^ may be different from the one in h e and 
h r . From the quark masses, which also contain e and e, one infers e ~ e 2 
[36]. 

From eq. (74) one obtains for the masses of charged leptons, light and 
heavy Majorana neutrinos, 



m e : rrin : m T ~ M\ : M2 : M3 



.5 . 



e:l, 



mi : m2 : m3 ~ e 15 : e 3 : 1 . 



(75) 
(76) 



Like in the example with SU (5) xU(1)f symmetry, the mass of the heaviest 
Majorana neutrino, 



M 3 



V 2 

m 3 



~ 10 15 GeV , 



(77) 



implies that B — L is broken at the unification scale Aqut- 

The f M — v T mixing angle is related to the mixing of the charged leptons 
of the second and third generation [36], 



sin 6 MT ~ y/e + e . 
This requires large flavour mixing, 

•(1>)^ 1/2 



A 



1 



(78) 



(79) 



In view of the unknown coefficients O(l) the corresponding mixing angle 
sinB^r ~ 0.7 is consistent with the interpretation of the atmospheric neu- 
trino anomaly as — v T oscillation. 

1 Note, that with respect to ref. [36], e and 1 have been interchanged. 
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It is very instructive to compare the two scenarios of lepton masses and 
mixings described above. In the first case, the large — v T mixing angle 
follows from a non-parallel flavour symmetry. The parameter e, which char- 
acterizes the flavour mixing, is small. In the second case, the large — v T 
mixing angle is a consequence of the large flavour mixing e. The U(l)p 
charges of all leptons are the same, i.e., one has a parallel family structure. 
Also the mass hierarchies, given in terms of e, are rather different. This 
illustrates that the separation into a flavour mixing parameter e and coef- 
ficients 0(1) is far from unique. It is therefore important to study other 
observables which depend on the lepton mass matrices. This includes lep- 
ton flavour changing processes and, in particular, the cosmological baryon 
asymmetry. 

5. Calculating the baryon asymmetry 

5.1. MAJORANA NEUTRINO DECAYS 

The heavy Majorana neutrinos, whose exchange may erase any lepton asym- 
metry, can also generate a lepton asymmetry by means of out-of-equilibrium 
decays. This lepton asymmetry is then partially transformed into a baryon 
asymmetry by sphaleron processes. The decay width of the heavy neutrino 
Ni reads at tree level, 

T Di = T(N l ^H 2 + l) + T {Ni - HI + l c ) = -L^t )^M; . (80) 

A necessary requirement for baryogenesis is the out-of-equilibrium condi- 
tion Tdi < H\t=Mi [42], where H is the Hubble parameter at temperature 
T. From the decay width (80) one then obtains an upper bound on an 
effective light neutrino mass [43, 44], 



T=M 1 



< 10~ 3 eV. (81) 

Here is the number of relativistic degrees of freedom, Mp = (SttGn)^ 1 ^ 2 — 
2.4 • 10 18 GeV is the Planck mass, and we have assumed g* ~ 100, t>2 ~ 
174 GeV. More direct bounds on the light neutrino masses depend on the 
structure of the Dirac neutrino mass matrix. If the bound (81) is satisfied, 
the heavy neutrinos Ni are not able to follow the rapid change of the equi- 
librium distribution once the temperature of the universe drops below the 
mass M\. Hence, the deviation from thermal equilibrium consists in a too 
large number density of heavy neutrinos N\ as compared to the equilibrium 
density (cf. section 2.1). 
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Figure 7. Tree level and one-loop diagrams contributing to heavy neutrino decays. 



Eventually, however, the neutrinos will decay, and a lepton asymmetry 
is generated due to the CP asymmetry which comes about through interfer- 
ence between the tree-level amplitude and the one-loop diagrams shown in 
fig. 7. In a basis, where the right-handed neutrino mass matrix M = h r (R) 
is diagonal, one obtains 



£l 



T(N 1 ^lH 2 )-T(N 1 ^l c Hl 
r{N 1 ^lH 2 )+T{N 1 ^l c H%] 

1 ^ X, ,\2 



1 

8^ 



E Im (huhl) 



M„ 



M? 



f\Tk +9 T75 



Ml 



(h v hVj i =2 ,3 

here / is the contribution from the one-loop vertex correction, 

'1 + x 



Mf 



;(82) 



f(x) = y/x 



1- (1 + x) In 



X 



(83) 



and g denotes the contribution from the one-loop self energy [45, 46, 47], 
which can be reliably calculated in perturbation theory for sufficiently large 
mass splittings, i.e., |Mj — M\\ 3> |Tj — 



g(x) = 



\fx 
1-x 



For Mi < M 2 , M 3 one obtains 
3 1 



16vr 



(jU) u 



E Im [(^t) 



i=2,3 



Ml 

M 



(84) 



(85) 



In the case of mass differences of order the decay widths one expects an 
enhancement from the self-energy contribution [48]. 

The CP asymmetry (82) leads to a lepton asymmetry [42], 



Y L = nL ~^ = k fl 
s g* 



(86) 
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Here the factor k < 1 represents the effect of washout processes. In order to 
determine k one has to solve the full Boltzmann equations [18, 19]. In the 
examples discussed in sections 4.1 and 4.2 one obtains n ~ 10 -1 . . . 10 -3 . 
Important processes are the AL = 2 lepton Higgs scatterings mediated 
by heavy neutrinos (cf. fig. 5) since cancellations between on-shell contri- 
butions to these scatterings and contributions from neutrino decays and 
inverse decays ensure that no asymmetry is generated in thermal equilib- 
rium [17]. Further, due to the large top-quark Yukawa coupling one has to 
take into account neutrino top-quark scatterings mediated by Higgs bosons 
[18, 19]. These processes are of crucial importance for leptogenesis, since 
they can create a thermal population of heavy neutrinos at high tempera- 
tures T > Mi. Alternatively, the density of heavy neutrinos may be gen- 
erated by inflaton decays [49] . Obviously, the requested baryon asymmetry 
can only be generated if the heavy neutrinos are sufficiently numerous be- 
fore decaying. 

Leptogenesis has been considered for various extensions of the standard 
model (cf. [16]). In particular, it is intriguing that in the simple case of 
hierarchical heavy neutrino masses the observed value of the baryon asym- 
metry is obtained without any fine tuning of parameters if B — L is broken 
at the unification scale, Agut ~ 10 16 GeV [50]. The corresponding light 
neutrino masses are very small, i.e., m V2 ~ 3 • 10~ 3 eV, as preferred by the 
MSW explanation of the solar neutrino deficit, and m vz ~ 0.1 eV. Such 
small neutrino masses are also consistent with the atmospheric neutrino 
anomaly [32], which implies a small mass m„ 3 in the case of hierarchical 
neutrino masses. This fact gave rise to a renewed interest in leptogenesis in 
recent years. 

5.2. BARYON ASYMMETRY AND BARYOGENESIS TEMPERATURE 

We can now evaluate the baryon asymmetry for the two patterns of neu- 
trino mass matrices discussed in sections 4.1 and 4.2. Since for the Yukawa 
couplings only the powers in e are known, we will also obtain the CP asym- 
metries and the corresponding baryon asymmetries to leading order in e, 
i.e., up to unknown factors 0(1). Note, that for models with a U(1)f gen- 
eration symmetry the baryon asymmetry is 'quantized', i.e., changing the 
U(1)f charges will change the baryon asymmetry by powers of e [37]. 

5.2.1. SU{5) x U(1) F 

In this case one obtains from eqs. (82) and (70), 



° 4 

Ei ~ e 

16vr 



(87) 
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z = M,/T 

Figure 8. Time evolution of the neutrino number density and the lepton asymmetry for 
the SU(5) x U(1)f model. The solid line shows the solution of the Boltzmann equation for 
the right-handed neutrinos, while the corresponding equilibrium distribution is represented 
by the dashed line. The absolute value of the lepton asymmetry Yl is given by the dotted 
line and the hatched area shows the lepton asymmetry corresponding to the observed 
baryon asymmetry. 

From eq. (86), e 2 ~ 1/300 (69) and g* ~ 100 one then obtains the baryon 
asymmetry, 

Y b ~k 10~ 8 . (88) 

For k ~ 0.1... 0.01 this is indeed the correct order of magnitude! The 
baryogenesis temperature is given by the mass of the lightest of the heavy 
Majorana neutrinos, 

Tb ~ Mi ~ e 4 M 3 ~ 10 10 GeV . (89) 

For this model, where the CP asymmetry is determined by the mass hierar- 
chy of light and heavy Majorana neutrinos, baryogenesis has been studied 
in detail in [50]. The generated baryon asymmetry does not depend on the 
flavour mixing of the light neutrinos, in particular the v a — v T mixing angle. 

The solution of the full Boltzmann equations is shown in fig. 8 for the 
non-supersymmetric case; the supersymmetric model has been studied in 
[51]. The initial condition at a temperature T ~ 10Mi is chosen to be a state 
without heavy neutrinos. The Yukawa interactions are sufficient to bring 
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z = M,/T 

Figure 9. Solution of the Boltzmann equations for the SU(3) c x SU(3)l* SU(3)rxU(1)f 
model. 

the heavy neutrinos into thermal equilibrium. The approach of the heavy 
neutrino number density to the equilibrium density and the evolution of the 
lepton asymmetry are analogous to GUT baryogenesis with heavy bosons 
[52]. At temperatures T ~ M\ the familiar out-of-equilibrium decays sets 
in, which leads to a non-vanishing baryon asymmetry. The final asymmetry 
agrees with the estimate (88) for k ~ 0.1. The dip in fig. 8 is due to a change 
of sign in the lepton asymmetry at T ~ M\. 



5.2.2. SU(3) C x SU(3) L x SU(S)r x U(1) f 

In this model the neutrino Yukawa couplings (74) yield the CP asymmetry 

which correspond to the baryon asymmetry (cf. (86)) 

Y B ~ k 10" 6 . (91) 

Due to the large value of e the CP asymmetry is two orders of magnitude 
larger than in the SU(<j) x U(1)f model. However, washout processes are 
now also stronger. The solution of the Boltzmann equations is shown in 
fig. 9. The final asymmetry is again Y# ~ 10~ 9 which corresponds to k ~ 
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1(T 3 . The baryo genesis temperature is considerably larger than in the first 

case, 

Tb ~ Mi ~ e 5 M 3 ~ 10 12 GeV . (92) 

The baryon asymmetry is largely determined by the parameter rh\ de- 
fined in eq. (81) [19]. In the first example, one has rh\ ~ m„. In the second 
case one finds fh\ ~ 7773. Since rn v and 777.3 are rather similar it is not too 
surprizing that the generated baryon asymmetry is about the same in both 
cases. 

6. Implications for dark matter 

The experimental evidence for small neutrino masses, the see-saw mecha- 
nism and the out-of-equilibrium condition for the decay of the heavy Ma- 
jorana neutrinos are all consistent and suggest rather large heavy neutrino 
masses and a correspondingly large baryogenesis temperature. For thermal 
leptogenesis models one typically finds [16], 

Tb ~ Mi > 10 7 GeV . (93) 

In the particularly attractive supersymmetric version of thermal lepto- 
genesis one then has to consider the following two issues: the consistency 
of the large baryogenesis temperature with the 'gravitino constraint' and 
the size of other possible contributions to the baryon asymmetry. A large 
asymmetry can in principle be generated by coherent oscillations of scalar 
fields which carry baryon and lepton number [8] . It appears likely, however, 
that the interactions of the right-handed neutrinos are sufficiently strong 
to erase such primordial baryon and lepton asymmetries before thermal 
leptogenesis takes place [53]. 

The couplings of gravitinos with matter are essentially model indepen- 
dent. Their cosmological effects therefore provide very interesting informa- 
tion about possible extensions of the standard model. It was realized long 




Figure 10. Typical gravitino production processes mediated by gluon exchange. 
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ago that standard cosmology requires gravitinos to be either very light, 
rriQ < IkeV [54], or very heavy, iuq > lOTeV [55]. These constraints are 
relaxed if the standard cosmology is extended to include an inflationary 
phase [56, 57, 58]. The cosmologically relevant gravitino abundance is then 
created in the reheating phase after inflation in which a reheating tem- 
perature Tr is reached. Gravitinos are dominantly produced by inelastic 
2^2 scattering processes of particles from the thermal bath. The grav- 
itino abundance is essentially linear in the reheating temperature Tr. It is 
intriguing that for temperatures Tr ~ 10 10 GeV, which are of interest for 
leptogenesis, gravitinos with mass of the electroweak scale, thq ~ 100 GeV, 
can be the dominant component of cold dark matter [59]. 

The bounds on the reheating temperature depend on the thermal grav- 
itino production rate which is dominated by two-body processes involving 
gluinos (g) (cf. fig. 10). On dimensional grounds the production rate has 
the form 

T(T) oc ± T 3 , (94) 

where M = (87tGat)~ 1/2 = 2.4 • 10 18 GeV is the Planck mass. Hence, the 
density of thermally produced gravitinos increases strongly with tempera- 
ture. 

The production rate depends on the ratio m g /mQ, the ratio of gluino 
and gravitino masses. The ten 2^2 gravitino production processes were 
considered in [57] for m g <C m,Q. The case m g » ra^, where the goldstino 
contribution dominates, was studied in [58]. Four of the ten production pro- 
cesses are logarithmically singular due to the exchange of massless gluons. 
The complete result for the logarithmically singular part of the produc- 
tion rate was obtained in [59]. The correct finite part can be obtained by 
means of a hard thermal loop resummation, which was first implemented 
in the case of axion production in a QED plasma [60] and recently also for 
gravitino production in a QCD plasma [61]. The result for the Boltmann 
collision term reads 




(N + n f ) + 0.5781n/| , (95) 



where 



m 2 g = ^(N + n f ) (96) 



is the thermal gluon mass squared; N is the number of colours and n/ is 
the number of colour triplet and anti-triplet chiral multiplets. The QCD 
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Figure 11. Contribution of gravitinos to the density parameter ft^h 2 for different grav- 
itino masses nig as function of the reheating temperature Tr. The gluino mass has been 
set to m~ g = 700 GeV [61]. 



coupling g(T) ~ 0.85 for T ~ 10 10 GeV. For the supersymmetric standard 
model with Nq = 3 and rif = 6 this implies m g > T. Hence, the usually- 
assumed separation of scales, g 2 T <C gT <C T, appears problematic and 
higher-order corrections may be important. 
Using the Boltzmann equation, 



(97) 



one can calculate the gravitino abundance at temperatures T < Tr, assum- 
ing constant entropy. One finds 



nrad{T) g*s(T R ) H(T R )n raA (T R ) ' 



(98) 



where g+s{T) is the number of effectively massless degrees of freedom 
[2]. For T < 1 MeV, i.e. after nucleosynthesis, g*s(T) = |f, whereas 
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9*s(Tr) = in the supersymmetric standard model. With H(T) = 
(^(T)7r 2 /90) 1 / 2 T 2 /M one obtains in the case of light gravitinos (m^ <C 
rrig(fi), (i ~ 100 GeV) from eqs. (98) and (95) for the gravitino abundance 
and for the contribution to S7/i 2 , 




Here we have used n rad (T) = C(3)T 3 /vr 2 , and m~ g (T) = g 2 (T)/g 2 (p)m g (p). 
The new result for ti^h 2 is smaller by a factor of 3 compared to the result 
given in [59]. This is due to a partial cancellation between the logarithmic 
term and the constant term in eq. (95). 

It is remarkable that reheating temperatures Tr = 10 8 — 10 10 GeV lead 
to values fl^h 2 = 0.01 ... 1 in an interesting gravitino mass range. This is 
illustrated in fig. 11 for a gluino mass m g = 700 GeV. As an example, for 
T R ~ 10 10 GeV, rriQ ~ 80 GeV and h ~ 0.65 [2] one finds ty* = 0.35, which 
agrees with recent measurements of the matter density $7m [2]. 

In general, to find a viable cosmological scenario one has to avoid two 
types of gravitino problems: For unstable gravitinos their decay products 
must not alter the observed abundances of light elements in the universe, 
which is referred to as the big bang nucleosynthesis (BBN) constraint. 
For stable gravitinos this condition has to be met by other super parti- 
cles, in particular the next-to-lightest super particle (NSP), which decays 
into gravitinos; further, their contribution to the energy density of the 
universe must not exceed the closure limit, i.e. Qq = Pq/p c < 1 5 where 
p c = 3H 2 M 2 = 1.05/i 2 10 -5 GeV cm -3 is the critical energy density. 

Consider first the constraint from the closure limit. The condition fig = 
YG m Q n xa,&l Pc < 1 yields a boundary in the m^-nig plane which is shown 
in fig. 12 for three different values of the reheating temperature Tr. The 
allowed regions are below the three solid lines, respectively. 

With respect to the BBN constraint, consider some nonrelativistic par- 
ticle X which decays into electromagnetically and strongly interacting rel- 
ativistic particles with a lifetime tx- Roughly speaking, the decay changes 
the abundances of light elements the more the longer the lifetime tx and 
the higher the energy density mxYx are. These constraints have been stud- 
ied in detail by several groups [62, 63, 64]. For most supergravity models 
they rule out the possibility of unstable gravitinos with run ~ 100 GeV for 
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Figure 12. Upper and lower bounds on the gluino mass and the NSP mass as func- 
tions of the gravitino mass. The full lines represent the upper bound on the gluino mass 
■nig > m N sp for different reheating temperatures from the closure limit constraint. The 
dashed line is the lower bound on itinsp which follows from the NSP lifetime [61]. 



Tr ~ 10 GeV, although even larger reheating temperatures are acceptable 
in some cases [65]. 

For stable gravitinos the NSP plays the role of the particle X. The 
lifetime of a fermion decaying into its scalar partner and a gravitino is 

m%M 2 

7nsp = 48tt^ . (101) 

m NSP 

For a sufficiently short lifetime, tnsp < 2 • 10 6 s, it is sufficient to require 
the energy density which becomes free in NSP decays to be smaller than 
m x Y x < 4 • 10~ 10 GeV, which corresponds to U x h 2 < 0.008. This con- 
straint can be satisfied since the NSP relic density is rather model depen- 
dent. For neutralinos in the MSSM the energy fraction Q x can vary over 
eight orders of magnitude (cf. fig. (13)). The lifetime constraint tnsp < 
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2 • 10 s yields a lower bound on super particle masses which is represented 
by the dashed line in the mQ-m N gp/g plane in fig. 12. 

In order to decide whether the second part of the BBN constraint, 
^nsp^ 2 < 0.008, is satisfied, one has to specify which particle is the NSP. 
The case of a higgsino-like neutralino as NSP has been studied in [59]. A 
detailed discussion for a scalar r-lepton as NSP has been given in [66], [67]. 

7. Outlook 

Recent developments in theoretical and experimental particle physics sup- 
port the idea of leptogenesis according to which the cosmological matter 
density has been created in decays of heavy Majorana neutrinos. On the 
theoretical side, detailed studies of the electroweak phase transition and 
sphaleron processes have shown that the matter-antimatter asymmetry is 
related to neutrino properties. On the experimental side, the solar and at- 
mospheric neutrino deficits have been observed, which can be interpreted 
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as a result of oscillations between three species of light Majorana neutrinos. 

It is very remarkable that these hints on the nature of lepton num- 
ber violation fit very well together with the leptogenesis mechanism. For 
hierarchical neutrino masses, with B — L broken at the unification scale 
Agut ~ 10 16 GeV, the observed baryon asymmetry ng/s~ 1CT 10 is natu- 
rally explained by the decay of heavy Majorana neutrinos. The correspond- 
ing baryogenesis temperature is Tg ~ 10 10 GeV. 

The consistency of this picture has implications in particle physics and 
cosmology. In unified theories the pattern of neutrino masses and mixings 
is related to lepton flavour and quark flavour changing processes. In su- 
persymmetric theories the mass spectrum of superparticles is constrained 
by the cosmological bound on the number density of gravitinos which may 
be the dominant component of dark matter. Further, the realization of the 
rather large baryogenesis temperature in models of inflation should have 
observable consequences for the anisotropy of the cosmic microwave back- 
ground. 
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